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Clusters of galaxies are the largest gravitationally bound objects in the Universe, con- 
taining about 10 15 solar masses of hot (10 8 K) gas, galaxies and dark matter in a typical 

P-i- 

volume of about lOMpc 3 . Magnetic fields and relativistic particles are mixed with the 
| ! g as as ^ h y g ia„t radi „ ha.oe, which arise from di^e, m „-sca, e sy , 
chrotron radiation at cluster center 1 2 . Radio haloes require that the emitting electrons 
are accelerated in situ (by turbulence) 3 6 , or are injected (as secondary particles) by 

oo : 

proton collisions into the intergalactic medium 7 10 . They are found only in a fraction of 
massive clusters that have complex dynamics 11 14 , which suggests a connection between 

oo : 

these mechanisms and cluster mergers. Here we report a radio halo at low frequencies 

> : 

associated with the merging cluster Abell521. This halo has an extremely steep radio 
spectrum, which implies a high frequency cut-off; this makes the halo difficult to detect 
with observations at 1.4GHz (the frequency at which all other known radio haloes have 
been best studied). The spectrum of the halo is inconsistent with a secondary origin of 
the relativistic electrons, but instead supports turbulent acceleration, which suggests that 
many radio haloes in the Universe should emit mainly at low frequencies. 

The turbulent re-acceleration model 5,6 ' 15-17 for the origin of giant radio haloes assumes 
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that fossil relativistic particles are re- accelerated by merger-induced turbulence to the energies 
necessary to produce the observed radio synchrotron emission in relatively weak magnetic 
fields. The acceleration of fast particles by turbulence, known to be an important process in 
astrophysics, is due to the resonant scattering of these particles by the turbulent waves (the 
energy spectrum of the turbulence can be thought of as the superposition of the contribu- 
tion from turbulent waves with different length scales), leading to a stochastic energization of 
particles and to the damping of the waves 18 . The connection between mergers and particle 
acceleration by turbulence in galaxy clusters is complex. It is argued 15,19 ' 20 that it should take 
place on a timescale of the order of a clustersubcluster crossing time (~lGyr), during which 
turbulence is continuously injected on scales of the order of the subcluster size, transported at 
smaller scales and then dissipated into heating of the intergalactic medium and acceleration 
of relativistic particles over a fairly large volume. This argument is supported by recent radio 
observations of a complete sample of X-ray luminous clusters that allow a clear separation 
between clusters with radio haloes and radio-quiet clusters, suggesting, from the fraction of 
clusters with radio haloes, that the particle acceleration mechanisms operate sporadically, on 
timescales <lGyr and in connection with cluster mergers 21 . Abell521 is an X-ray luminous 
(8.2 x 10 37 W in the 0.12.4-keV band) and massive (~ 2 x 10 15 solar masses) galaxy cluster at 
redshift z=0.247 with ongoing multiple merging episodes 22 ' 23 . Here we report the discovery 
of a giant radio halo in this cluster by means of deep observations with the Giant Metrewave 
Radio Telescope (GMRT, India) at 240, 325 and 610MHz. In Fig. la we show the radio 
image of Abell521 at 240MHz, where the radio halo is best imaged. To highlight the diffuse 
emission, point sources visible in the full-resolution images were subtracted in producing the 
low-resolution images in Fig. 1. The radio halo is coincident with the cluster X-ray emit- 
ting region and correlates with the X-ray emission of the hosting cluster, which is typical of 
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other radio haloes 1,2 (see also Supplementary Information). Previous higher frequency Very 
Large Array (VLA) observations at 1,400MHz did not find this diffuse emission, instead ob- 
serving only the radio relic located on the southeastern boundary of the cluster 23 . The relic 
coincides with a possible shock front, generated by recent in-fall of a subcluster along the 
northwest/southeast direction, where relativistic electrons are currently accelerated 24 . Figure 
1 clearly shows that the radio halo becomes progressively more dominant over the radio relic 
at lower frequencies, indicating that its spectrum is much steeper than that of the relic, which 
has a 1.5 (ref. 24; flux is proportional to v~ a , where v denotes frequency). The patchy 
structure of the radio halo at 610MHz (Fig. lb) indicates the observational difficulty in imag- 
ing the emission, the surface brightness of which is already fading at this frequency (note that 
the largest diffuse emission detectable with the GMRT at 610 MHz is ps 0.28 degrees, e.g. 
ref. 14). The halo disappears between 610 and 1,400MHz; only an upper limit on the flux of 
the radio halo at 1,400MHz can be derived, although faint residual emission in the cluster is 
still present at this frequency (Fig. lc and Supplementary Information). The flux densities of 
the radio halo at 240, 325 and 610MHz are plotted in Fig. 2 together with the upper limits 
at 74 and 1,400MHz. The upper limits were evaluated by injecting fake radio haloes with 
different flux densities into the observed data sets, following ref. 21, to estimate the sensitiv- 
ity of the observations to diffuse emission on the halo length scale. In particular, the upper 
limit at 74MHz was derived from VLA Low-frequency Sky Survey data 25 , whereas the upper 
limit at 1,400MHz was derived from the analysis of the archival VLA data (Supplementary 
Information). The important result is that the average value of the spectral index, a ps 2.1, 
is much larger than that of any other known radio halo (typical spectral index is a ps 1.2-1.3 
(refs 1, 2)). These extreme spectral properties make Abell521 a unique system for addressing 
the origin of the emitting particles in radio haloes. Such an extremely steep spectrum and 
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the downward spectral curvature (Fig. 2) imply a spectral cut-off at high frequency, which 
is a well-known signature of turbulent acceleration 4-6,16 ' 17,19 . Synchrotron theory implies a 
corresponding cut-off in the spectrum of the emitting electrons at E e « 1.45 nr 1 ^ 2 (v c /300) 1//2 
GeV, where v c is the cut-off frequency (measured in megahertz) and B n j is the magnetic field 
(measured in nT). In the framework of the turbulent acceleration scenario, E e pinpoints the 
energy at which the timescale of the electron radiative losses becomes equal to that of the tur- 
bulent acceleration. The timescale of electrons emitting at v c is estimated taking into account 
the red shift-dependent inverse-Compton losses against the cosmic microwave background and 
the synchrotron losses 26 , as follows, where v c is again measured in megahertz: 



B l J T 2 (v c /300)- 1 / 2 

-r ~ o 95 nT - Gvr (1) 

U - y '(l+ Z ) 4 + (5„ r /0.32)2 Uyr {l) 



This means that the electrons responsible for the observed emission should be accelerated 
on a timescale of ~ (1.1 — 1.4) x 10 8 years, for 0.1 - 0.5-nT magnetic fields in the radio halo 
region. Assuming that fast magnetosonic waves are responsible for the acceleration of the 
emitting particles, following ref. 20 we find that this acceleration efficiency can be achieved 
under the reasonable assumption that the energy density of these turbulent waves is 12-18% 
of the thermal energy (Fig. 2). 

Appreciable synchrotron emission can also be produced by secondary electrons injected 
by collisions between long-lived relativistic protons accumulated in the cluster and the thermal 
protons in the intergalactic medium, and secondary models have been proposed as alternatives 
to the re- acceleration model to explain radio haloes 7-10 . The very steep spectral slope of 
this radio halo rules out secondary models by means of a straightforward energy argument. To 
explain the spectrum of the radio halo through synchrotron radiation from secondary electrons, 
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the primary protons must have a very steep spectral energy distribution (N(p) <x p °, where 
<5 ps 4.2 and p denotes the particle momentum). The energy density of relativistic protons, £ p , 
required to match the observed synchrotron flux through a secondary model can be estimated 
following the formalism in ref. 27. For an average number density of thermal protons of 
vith ~ 1 > 500m -3 in the region of the radio halo (consistent with the average thermal density 
in the same region derived from X-ray observations 22 ) and the synchrotron flux measured 
at 325MHz, we find that £ p ranges from approximately 3 to 100 times the energy density 
of the thermal plasma for magnetic field values (averaged in the region of the radio halo) 
ranging from B=0.5nT to B=0.1nT, respectively. This yields only a lower limit on the energy 
density of high-energy protons, because for 8 > 3 an additional (even dominant) contribution 
to the energy comes from suprathermal particles with kinetic energies <lGeV. A secondary 
origin for the emitting electrons thus implies the unrealistic situation in which clusters are 
dominated by non-thermal protons. It also violates present upper limits on the energy density 
of these particles derived from y-ray observations of several clusters, at the 20% level 28 ' 29 . 
If the spectrum of these protons has 6 < 4, then the synchrotron signal from the secondary 
electrons produced by proton-proton collisions cannot exceed the 1.4-GHz upper limit in Fig. 
2, placing corresponding limits on the energy density of the primary protons. Figure 3 shows 
upper limits on the energy density of the primary protons in the region of the radio halo as a 
function of the magnetic field. These limits were obtained following ref. 21 and show that the 
energy density of relativistic protons in the cluster is less than 1% of the thermal component 
for B > 0.2nT and 6 < 2.5. A similar conclusion has recently been reached in the analysis of 
a statistical sample of galaxy clusters without radio haloes 21 . However, in that case a larger 
energy content of protons was still possible by assuming that clusters without radio haloes 
have magnetic fields much smaller than those with radio haloes. This alternative (ad hoc) 
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possibility can be reasonably ruled out in our case because Abell521 hosts a radio halo with 
bolometric radio luminosity comparable to that of classical radio haloes (vP(v) 10 34 W). 
Future observations with the Fermi Gamma-ray Space Telescope (formerly the Gamma-ray 
Large Area Space Telescope) will reveal galaxy clusters in y-rays in cases where the energy 
content of relativistic protons is significantly larger than about 1% of the thermal plasma. The 
combination of this future data and limits in the radio band (Fig. 3 and ref. 21) may thus 
provide a powerful tool for constraining the magnetic field strength in galaxy clusters. 

As we look at applying the turbulent re- acceleration model to other galaxy clusters, it 
should be stressed that the maximum energy to which electrons can be re- accelerated and, 
ultimately, the cut-off frequency in the spectra of radio haloes depend on the level of turbu- 
lence and on the properties of the turbulent waves. The spectral cut-off affects our ability to 
detect radio haloes in the Universe, introducing a strong bias against observing them at fre- 
quencies substantially larger than v c . Currently known radio haloes are observed at gigahertz 
frequencies, requiring an efficient turbulent acceleration mechanism. These haloes must re- 
sult from the rare, most energetic merging events and are therefore hosted only in the most 
massive, hottest clusters 19,30 , in line with observations u ~ 14 . On the other hand, the majority 
of radio haloes should form during much more common, but less energetic, merging events, 
for example between a massive cluster and a substantially smaller subcluster (with mass ra- 
tio >5) or between two similar clusters with mass < 10 15 solar masses J '. However, these 
sources, with a cut-off in the synchrotron spectrum at v c < 1GHz, should be visible only at 
lower frequencies, because their spectrum should be similar to that of the low-frequency radio 
halo in Abell521. Future high-sensitivity radio telescopes operating at low frequencies, such 
as the Low Frequency Array (LOFAR) and the Long Wavelength Array (LWA), are expected 
to discover the majority of these sources and also to test their connection with cluster mergers. 
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Figure 1: Radio and X-ray images of Abell 521. Low resolution radio contours (after subtraction 
of discrete sources in the field identified at full resolution) overlaid on the Chandra 0.5-4 keV X-ray 
image (corrected for background and exposure and smoothed with a a=6 arcsec Gaussian, as in ref. 24). 
Panel a): GMRT 240 MHz contours at 35 arcsec x 35 arcsec resolution. The r.m.s. noise level in the 
cluster region is 220 /iJy/beam (full resolution used for subtraction of discrete sources has 15.6 arcsec 
x 12.3 arcsec beam and r.m.s. noise =190 /iJy/beam). Contours are at -0.66 (red), 0.66, 1.3, 2.6, 5.2, 
7, 10, 14, 19 mJy/beam. Panel b): GMRT 610 MHz contours at 35 arcsec x 35 arcsec resolution. The 
r.m.s. noise level in the cluster region is 58 iiJy/beam (full resolution used for subtraction of discrete 
sources has 9 arcsec x 4 arcsec beam and r.m.s. noise =35 iiJy/beam). Contours are at -0.17 (red), 0.17, 
0.34, 0.48, 0.67, 0.94, 1.3, 1.9, 2.6, 3.7, 5.2 mJy/beam. The dashed circle indicates the region in which 
fluxes of the radio halo at different frequencies (reported in Fig. 2) are measured. Panel c): VLA 1.4 
GHz contours at 25 arcsec x 25 arcsec resolution. Here only the discrete sources within 3 arcmin from 
the cluster centre are subtracted (other images are presented in Supplementary Information together 
with a discussion on the subtraction of discrete sources at this frequency). The r.m.s. noise level in 
the cluster region is 26 /iJy/beam (the full resolution used for subtraction of discrete sources has 12.7 
arcsec x 6.9 arcsec beam and r.m.s. noise =15 iiJy/beam). Contours are at -0.08 (red), 0.08, 0.16, 0.22, 
0.32, 0.45, 0.63, 0.89, 1.25, 1.76, 2.48, 3.5, 4.9 mJy/beam. Coordinate system, J2000. 

At the moment this connection cannot be tested, owing to the lack of observations of samples 
of galaxy clusters at low radio frequencies. 
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Figure 2: The spectrum of the radio halo. Data are overlaid on secondary (long-dashed line) and 
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Figure 3: Limits to the energy density of relativistic protons. Curves give the upper limit to the 
ratio between the energy density of relativistic protons and that of thermal gas, e t h, as a function of the 
magnetic field in Abell 521. Calculations are shown assuming = 1,500 m~ 3 and ksT =7 keV. From 
bottom to top we adopt a spectra of the primary protons : 8= 2.3, 2.5, 2.7, 2.9 (N(p) oc p~ s ). 
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To complement the information given in the 
main body of the Letter, these Sections pro- 
vide supplementary discussion on the methods 
adopted to derive upper limits to the flux density 
of the halo reported in Figure 2 (Sect. 1) and on 
the origin of the halo (Sect. 2 and 3). 

Upper limits to the flux of the radio halo 
at 1.4 GHz 

The diffuse emission in the center of Abell 521 dis- 
covered by the low frequency radio observations pre- 
sented in this Letter was not detected in VLA obser- 
vations at 1.4 GHz! 17 ]. We re-analyzed those BnC 
VLA data taken from the archive. Given the sour- 
thern declination of the cluster, they can be consid- 
ered equivalent to a C configuration observation. The 
quality of the data is generally good, except for oc- 
casional intermittent RFI in IF2 (1.44 GHz) and at 
the very end of the observation (both 25-MHz IFs 
affected), when the source and its secondary cali- 
brator were observed at low elevation and data on 
a few short baselines had corrupted amplitudes (very 
high circular polarization). RFI excision was carried 
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out with great care. We also paid particular atten- 
tion to the analysis of the secondary calibrator data, 
whose visibility turned out to be affected by substan- 
tial contribution from sources in the field. Therefore 
we first followed the standard calibration procedure, 
then we used the image of the secondary calibra- 
tor obtained in this way to improve the a-priori am- 
plitude (and phase) calibration accuracy. The final 
dataset was obtained after a number of phase-only 
self-calibrations. Solutions had generally very high 
signal to noise ratio. 

The final images were obtained by considering a 
central square, 0.8° wide, with additional 10 smaller 
(0.1°) facets on confusing sources. In the full reso- 
lution image (12.7 arcsec x6.9 arcsec in p. a. 75°) 
the r.m.s. noise level is 15 /xJy/beam. We also pro- 
duced images with lower resolution by applying nat- 
ural weighting to the data and also using a restricted 
uv-range. 

A low resolution 30 arcsec x 30 arcsec image of 
the region of Abell 521 is shown in Figure SI. Sus- 
picious, very faint, diffuse emission appears in the 
region of the radio halo discovered at lower frequen- 
cies; however, most of this emission is due to the 
convolution of point sources in the same region (the 
image after subtraction of point sources is given in 
Figure 1) and the detection of the radio halo in Abell 
521 at 1.4 GHz would not have been firmly estab- 
lished from these observations only. 

We thus estimated the upper limit to the flux den- 
sity of the radio halo in these data at 1.4 GHz. This 
should account also for the possibility that a fraction 
of the flux density on the halo scale may be lost due 
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Figure 4: Low resolution images of Abell 521 at 1.4 GHz. The upper panel shows the image from 
our analysis of the original data, while lower panels show images from the analysis of the data containing 
injected fake radio halos with ©50% =2.15 arcmin and F= 4 (left) and 6 (right) mjy. Crosses mark the 
position of the point sources identified in the high resolution image in the central region of Abell 521 and in 
the region where fake radio halos have been injected. Radio contours are at 80, 160, 240, .. jUJy/beam level; 
negative contours at -80/xJy/beam level are reported in dashed lines. The r.m.s. noise level in the images is 
in the range 28-30 jUJy/beam and the restoring beam is 30 arcsec x 30 arcsec. 
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to the uv coverage at short baselines of the obser- 
vations. In order to derive a solid and conservative 
upper limit we follow! 715 !. We injected fake halos 
in the uv-components of the dataset by means of the 
MPS task UVMODl 31 ! and imaged the datasets with 
the same set of parameters applied to the original 
data. This procedure allows to evaluate the response 
of the 1.4 GHz VLA observation to the presence of 
diffuse emission on the halo scale, as the uv com- 
ponents of the fake radio emission are sampled by 
means of the real uv coverage of the observation. 

In particular, we model the brightness profile of 
the fake halos with sets of optically thin spheres with 
different radius and flux density, and obtain families 
of fake radio halos with total flux densities ranging 
from 3 to 15 mJy, with sizescale enclosing 50% of 
the flux @ 50 % =1.9-2.3 arcmin (« 440 and 530 kpc, 
respectively) and maximum angular extension « 5.2 
arcmin (ps 1.2 Mpc). 

Figure S 1 shows the results obtained for fake ra- 
dio halos with ©50% =2.15 arcmin and injected flux 
density F= 4 and 6 mJy. In the case of the halo with 
F= 4 mJy, about half of the injected flux is recovered 
by the analysis of the dataset, and the resulting dif- 
fuse emission in the image is still brighter than that 
at the center of Abell 521. 

The analysis of the families of datasets with fake 
radio halos allowed us to derive a conservative up- 
per limit to the flux density at 1.4 GHz of the radio 
halo in Abell 521, F < 5 mJy. This should be com- 
pared with the residual diffuse emission at F ps 1 — 2 
mJy level measured in the region of the radio halo af- 
ter subtraction of discrete sources (Figure SI) in our 
images. Because the flux density of these sources, 
F = 3 — 4 mJy, is not negligible at 1 .4 GHz, they 
were fitted in the full resolution image (with beam 
= 12.7 arcsec x 6.9 arcsec and r.m.s. =15 jUJy/beam) 
by using the AIPS task JMFIT solving for a Gaussian 
component (which has then been subtracted out) and 
a non-zero level (i.e. diffuse emission, which is not 
subtracted). 

Projection effects and the radio halo 

Ghost radio plasma, confined within the intergalactic 
medium, may be revived by the passage of a shock 
wave through adiabatic compression 32 ; indeed this 
mechanism has been proposed to explain radio relics 



as a scenario alternative to shock acceleration. 

In principle, a large scale shock crossing Abell 
521 along the line of sight may compress ghost radio 
plasma producing detectable synchrotron radiation 
with steep spectrum projected on the cluster central 
region, that may be classified as radio halo. How- 
ever, the steep spectrum of the halo in Abell 521, its 
large {k, Mpc) extension and - in particular - the mor- 
phological correlation between radio halo and the X- 
ray emission from the cluster can be used to conclude 
that this scenario is very unlikely. 

In this context we have to assume that the bub- 
ble of ghost plasma extends on the Mpc scale. This 
implies a dynamical age of this plasma > Gyr dur- 
ing which relativistic electrons lose energy via ra- 
diative and adiabatic losses. This results in a cut- 
off in the synchrotron spectrum emitted by the ghost 
plasma at very low frequencies, Vb,o <MHz, before 
the shock-compression phase. The passage of the 
shock compresses the ghost plasma by a factor 32 
x = {Pds/Pus? 1 *, where P ds and Pus are the pres- 
sures measured in the thermal medium around the 
ghost plasma downstream and upstream the shock, 
respectively. The adiabatic compression boosts the 
frequency of the cut-off of the synchrotron spec- 
trum at 33 v c v c o x 4 / 3 ~ Pds/Pus- In order to have 
V c ~ 200 — 300 MHz as measured in the spectrum 
of the radio halo in Abell 521, this would require 
Pds/Pus > 300 implying a Mach number of the shock 
(from the relation between pressure jump and Mach 
number 34 ) : 

M ^( 4 tH' /2>15 {SA) 

Strong large-scale shocks with M > 10 are rare 
in galaxy clusters 35 - 36 - 37 and may occur only outside 
the virialised region of the cluster where the temper- 
ature of the gas is well below 1 keV. Thus the com- 
pression of the ghost plasma should happen at > 3 
Mpc distance from the center of Abell 521 (provided 
that a ghost-Gyr old plasma at such distances man- 
ages to survive from mixing with the intergalactic 
medium). 

At the same time the radio halo is correlated with 
the cluster X-ray emission and its radio brightness 
decreases by a factor «6-8 on the scale of the core 
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of Abell 521. In order to reproduce this morphologi- 
cal connection by chance the ghost radio bubble, at a 
distance of > 3 Mpc from the cluster center, should 
be aligned with the center of Abell 521 within a pro- 
jected distance of a cluster core-radius, m 300 kpc. 
Geometrically, this implies that the direction con- 
necting Abell 521 and the ghost bubble should make 
an angle B < 1° to the line of sight and the chance 
probability for that, « 1 — cos(0), is < 1%. 

The ad hoc assumptions and this low probability 
allow us to conclude that this scenario is very un- 
likely. 



Gyr (Eq.l), and thus we neglect adiabatic losses of 
electrons. This last assumption is justified since 
the gas in the downstream region, over-pressured 
with respect to the ambient gas, should expand in 
a timescale that is a fraction of m <p r /c s 0.5 Gyr, 
where r ps 500 — 700 kpc is the longitudinal size of 
the relic 18 and c s the sound velocity. 

Under these assumptions, the size of the tail at 
a given frequency can be estimated as l ta u k, v^T r , 
where x r is the radiative lifetime of electrons emit- 
ting at that frequency, and v^ s is the downstream ve- 
locity of the gas (in the shock frame, from the jump 
conditions 34 ) : 



On the size of the radio relic in Abell 521 

The radio halo in Abell 521 appears connected with 
the radio relic (Figure 1). This is similar to that ob- 
served in the case of Abell 520 33 and Abell 2256 38 ; 
a bridge connecting the halo and the relic is also ob- 
served in the Coma cluster 39 . 

Radio relics are interpreted as the result of the 
acceleration of relativistic particles due to large scale 
shocks that form in the intergalactic medium during 
cluster mergers or accretion of matter 40 41 . The relic 
in Abell 521 is found to coincide with a possible 
shock front with Mach number M « 2 — 4, generated 
by the recent infall of a subcluster along the North- 
West/South-East direction! 18 ] . 

Electrons accelerated at the shock flow away into 
the intergalactic medium in a tail of emitting plasma, 
and one might expect that part of the emission in the 
region of the radio halo is due to this tail. However 
in the following we show that this is unplausible : 
the extension of the tail at a given frequency is con- 
strained by the lifetime of the electrons emitting at 
that frequency (that depends on inverse Compton and 
synchrotron losses, and on adiabatic losses), which 
turns out to be much smaller than the >Mpc pro- 
jected distance between the relic and halo region. 

We adopt the most favorable geometry assum- 
ing that the shock is moving on the plane of the sky 
along the North- West/South-East direction, in which 
case the distance that separates the relic and halo re- 
gion is not affected by projection effects. We also as- 
sume that expansion in the shocked gas plays a role 
on a timescale longer than the radiative timescale 
of the electrons emitting at 240 MHz, that is 0.2 



where c" is the sound velocity in the intergalac- 
tic medium upstream of the shock. For a given Mach 
number of the shock c" can be estimated from the 
temperature of Abell 521 (taken as downstream tem- 
perature, and from the relations between temperature 
jump and Mach number 34 ) : 

r~ 16M 2 

C *~ 170 °\ W-l)(^ + 3) (to/s) <"> 

Combining S.2-S.3 and adopting z r « 0.2 Gyr 
for the electrons emitting at 240 MHz (Eq. 1), we 
obtain l ta u ps 200 kpc if the relic in Abell 521 in gen- 
erated by a shock with M rs 2 — 4. 

Remarkably, l ta u « 200 kpc agrees with the trans- 
verse size of the relic at 240 and 330 MHz, and with 
the size of the synchrotron spectral steepening in the 
relic as measured from radio observations 18 . 
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